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ABSTRACT 

The dynamics of Josephson-like oscillations between two coupled Bose-Einstein condensates is 
studied using the time-dependent variational method. We suppose that the quantum state of 
the condensates is a gaussian wave-packet which can translate and perform breathing shape 
oscillations. Under this hypotheses we study the influence of these degrees of freedom on the 
tunneling dynamics by comparing the full-model with one where these degrees of freedom are 
"frozen" at its equilibrium values. The result of our calculation shows that when the traps are 
not displaced the two models agree, whereas when they are, the models differ considerably, the 
former being now closer to its linear approximation. 
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The investigation of collective excitations of Josephson-coupled two-component Bose- 
Einstein condensates has been the subject of many papers 0-0. The main interest in those 
papers is the study of the dynamics associated to the exchange of atoms between the two con- 
densates. The basic hypotheses in most of these studies is that this exchange is "coherent" [|J, 
that is, without changing the quantum state of each condensate, which leads to a dynamics 
involving only the degree of freedom associated to the relative population and phase of the 
condensates. Based on this "frozen" model, a host of interesting phenomena is predicted, for 
example, the existence of a "phase-transition" to a non-symmetrical equilibrium composition 
in the limit of strong repulsion between the two condensates and "weak" Josephson coupling 




One exception, however, is the work of reference || where the effect of the change of the 
quantum state is taken into account. They use standard mean-field theory to derive coupled 
time-dependent Gross-Pitaevskii equations for the time-evolution of the order parameter of each 
condensate. They simplify the problem by treating the system in only one spatial dimension and 
they observe the effect of the mean-field on the non-linear oscillations of the relative population 
between the two condensates. 

In this paper we study the effect of the change of the quantum state in the tunneling 
dynamics. Our approach is in the spirit of reference |§, however we take a more qualitative 
point of view along the line of references We suppose that the quantum state of the system 

is a gaussian wave-packet which can translate and perform breathing shape oscillations. Under 
this hypotheses, we study the influence of these degrees of freedom on the tunneling dynamics 
by comparing the full-model with the limit where these degrees of freedom are "frozen" at its 
equilibrium values. In our study, we also investigate, in the context of the full-model, under 
what conditions the linear dynamics is valid. 

Our starting point is the action 
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S = J d 3 fdt [ ityj t)ipj (f, t) - S(r, t) 



(1) 



where i/jj(r,t), j = 1,2 are the order parameters of each condensate and the energy density is 
given by 



£(r,i) = $>;(r,0* 
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In the above expression V t J rap (r) is the trapping potential of each component which we take as 



V 3 

" trap 



m r 



uj 2 ± (x 2 + y 2 )+uj 2 z (z + lj z ) 2 



(3) 



where 71 = 1 and 72 = —1, Xjk are proportional to the scattering lengths ctjk for the collisions 
of atoms j and k of equal mass, Xjk = ^&jk, the term that depends on S takes into account the 
effects of detuning, and the last term is the Josephson interaction term where Q is the tunneling 
intensity ||. 

The condition that the action is stationary with respect to variations of the order parameters 
leads to the coupled mean-field equations for ipj(r, t) |§. In this paper we follow references |5|-|8[ 
and instead we parametrize the time dependence of the order parameter through the variation 
of a set of 2n parameters 0] which we denote by X = {Xi, X 2 , X 2n }, 



^•(r-,t)=^(f,X(t)) 



(4) 



Given this parametrization the action reduces to a "classical action" in terms Xj(t) 



where 



s = Jdt [£rf (x)l fc -£(x) 



(5) 



(6) 
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In (6) we introduced the state vector \ipj(X(t))) whose norm is equal to the population of the 
component j of the condensate and whose wave function is given by equation (4), and -E'(X) is 
the spatial integral of the energy density 

E{X) = [ d 3 fS(r,t) (7) 



with the order parameter parametrized as in eq. (4). 

Requiring that the action is stationary with respect to variations of the parameters leads to 
Hamiltonians-type equations of motion for Xi(t): 

E r $**=|^. (8) 
where the antisymmetric matrix is given by 

r (2) = _ (Q) 

jk dXj dX k 1 ; 



To proceed we should specify our choice of parameters. Previous work have shown 
that a parametrization which reproduces very well the equilibrium properties and collective 
excitations of a single- component condensate is one where the condensate density is a gaussian 
wave-packet which can translate and perform quadrupole and monopole shape oscillations and 
its corresponding motions. The study of the two-component condensate dynamics with the 
above parametrization is feasible and is under investigation. Here, in the spirit of reference |9, 
we simplify the problem by treating the system in only one spatial dimension. 

Thus we write the order parameter as: 

\^) = Jmje-^\ ( j )3 ) (10) 

where we put explicitly the degree of freedom associated to the population and phase of each 
condensate. The condensate orbitals \4>i) (1(J are taken as gaussian wave-packets which can 
translate and perform breathing shape oscillations, 



fafaX)) = (-=-) e > K) (11) 

where the parameters (Qj,Pj) are related to the translational degrees of freedom and (qj,Pj) 
to the "breathing" oscillation. Notice that in this parametrization the quantum state of each 
condensate (called here the condensate orbital) \4>j) does not depend explicitly on the phase 
and population of the condensates. 

Thus, we have twelve parameters which leads to a twelve dimensional "phase space" . Since 
E(K) depends only on the difference of phase between the two condensates, the equations of 
motion (8) give that the total number of particles is conserved. This reduces the number of 
degrees of freedom to ten, which we take as besides the g,, Qi,Pi, Pi, i = 1,2, the relative phase 
and relative population fraction. Thus, in the full-dynamics, we solve ten coupled first order 
equations, given an initial configuration. 

Since detr^ 2 ^(X) ^ 0, the equilibrium configuration of the system is determined by the 
condition 

dE 



dX j 



(X ) = 0. (12) 



Linearizing the equation (8) in the neighborhood of the equilibrium configuration, one has 

?^W*'W = ?a0ti^ (13) 

where x is the displacement from equilibrium 

x k = X k - X k0 (14) 

and in these equations denote that the quantities are evaluated at the equilibrium configu- 
ration. The solution of the coupled linear equations (13) define the small oscillation dynamics 
around this point. 

Finally, the frozen model is obtained by constraining all the degrees of freedom, except the 
relative population fraction and the relative phase at its equilibrium values. When this is done, 
the equations of motion (7) reduces to [|1],0 



2 ^(1-^)1/2 

r] = uj r {1 - r} 2 ) 1/2 sin^ (16) 

where 6 = 9 2 — Q\ and r\ = r] 2 — T]i, and in equations (15) the time is measured in unit of u , 
with u the trap frequency. The parameters defining the frozen model are : 

A = e 2 -e 1 -5 (17) 
A = An + A 22 - 2A 12 (18) 

90 r 

UJ R = T— ^(f, 0)Mr, 0) (19) 

(20) 

where e J is essentially the energy per-particle of each condensate (in units of ujq) 

N 



e 3 



\t + VLMi) + -j^j J d'fl^f, 0)| 4 ] /hw (21) 



and the Xjk are 



Xjk = N^- J d 3 r\Ur,0)\ 2 \Mr,0)\ 2 (22) 



We apply the formalism presented above to the system considered in reference 0, trapped 
87 Rb atoms in hyperfine states |/ = l,m/ = — 1) and |/ = 2,m/ = 1). Following this reference, 
the initial configuration is such that the only degree of freedom which is not at equilibrium is 
the phase difference and the detuning is chosen in such a way that the system is initially driven 
resonantly 0, which implies that A = 0. 

As a first example we will discuss the case where the traps are not displaced (zq = 0) and 
On = a 22 = 0,12 = ^Rb- In this case we have a complete agreement between the "frozen" and 
"full model", indicating that the tunneling dynamics is such that the quantum state of the 
condensates do not change. This result is a consequence of the fact that, in this symmetrical 
case, qi = q i0 , Qi = Q i0 , p i: = Pi = 0, i = 1, 2 is an invariant subspace (also called maximally 



decoupled subspace) of the system, that is, if the system is in this surface, it remains there at 
all time 0. 

We can also examine the importance of non-linear effects in this symmetrical case. The 
linearization of the frozen model gives a frequency of oscillation equal to the Rabi frequency 
ur, that in our case is equal to ujr = 0.100, which coincides with the frequency of the lowest 
energy normal mode of equations (13), u;^ = 0.100. When compared to the full model, the 
linear approximation basically differs only in the amplitude of the oscillation, being larger for 
the latter. 

As a more realistic application, we consider the case where the trap is not displaced (zq = 0) 
and the scattering lengths are equal to the experimental values as in Table 1 of reference || . In 
fig. 1 we show the results of the calculation corresponding to the frozen model, full-model and 
its linearization, when 0(0) = 0.77T as explained in the caption. As shown in these figures, there 
is again complete agreement between the "full" and "frozen" model. The translational degrees 
of freedom remain at its equilibrium values, whereas the displacements of the ones associated 
to the "breathing" oscillation are negligible (notice the scale of the figs. Id and le). Therefore, 
we can conclude that the invariant subspace found in the symmetrical case survives and, as a 
consequence, we have a complete agreement between the two models. 

In the case under investigation, the parameters of the frozen model are ojr = 0.100 and 
A = —0.017, showing that we are in the "strong Josephson" coupling regime > 1 fljj,|l2|- 
The linearization of the frozen model gives a frequency equal to up = ^ojr{ujr + ^) which 
coincides with the lowest energy normal mode of equation (13), whose frequency is = 0.096. 
As in the symmetrical case, the linear approximation of the full model differs from the other 
two models in the amplitude of the oscillation, with the former having a larger amplitude. Of 
course, as #(0) approaches the equilibrium value, #(0) = tt, this difference is less pronounced. 

Now we will turn to the investigation of the effect of the displacement of the traps on the 
system. To this end, we repeat the previous calculation with z = 0.15^0 JjJ. In fig. 2 we 
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show the curves for the relative population fraction, translation and "breathing" oscillation of 
each condensate component in the three models, with 0(0) = 0.77r as explained in the captions. 

From these figures we see that there is a good agreement between the full-model and its 
linear approximation whereas they both differ from the "frozen" model. From figure (2a) we 
observe that, for the relative population, the "frozen" model gives a much smaller amplitude and 
a much higher frequency, compared to the other two models. One the other hand, the results of 
the calculation according to the full model and its linear approximation are very close, the latter 
having a slighter higher amplitude and frequency. We observe the same qualitative pattern for 
the translational degrees of freedom, which oscillate in phase and the "breathing" ones, which 
oscillate out of phase. 

In the present case, the parameters of the frozen model are A = 128.4 and ujr = 0.0523, 
showing that we are in the regime of weak Josephson coupling -ffi < 1. The decrease of ujr and 
the increase of A can be understood as basically an effect of the displacement of the traps. This 
displacement diminishes the overlap between the condensate orbitals <f>j(r, 0), which decreases 
ujr and A 12 and increases A. 

The linear approximation of the frozen model gives a frequency of up ~ 1.832 much higher 
than the observed oscillation of the relative population. On the other hand, the lowest energy 
normal mode of the linear equations (13) has a frequency equal to 0.4466, which compares 
well with the frequency of oscillation for the relative population and the translation degree of 
freedom calculated in the linear approximation to the full-model. This result indicates that 
these degree of freedom has components basically in this normal mode. In turn, the time- 
evolution of the "breathing" variables, figs. (2d) and (2e), although dominated by this mode, 
show the presence of higher frequency modes. 

Comparing the two calculations at different displacements, we see that the full model pre- 
dicts that, when we displace the trap, the amplitude diminishes and the frequency increases, in 
agreement with reference ||. The frozen model gives the same qualitative prediction, however 
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it overestimates the size of the effect. 

To summarize, we have investigated the influence of the change of the quantum state in the 
tunneling dynamics of Josephson coupled two-component BEC, for different trap geometries. 
From our calculations, we have observed that when the trap is not displaced the tunneling 
dynamics is "coherent" , that is, the orbitals of each condensate do not change. As a consequence, 
we have agreement between the frozen model and the full model. On the other hand, when 
we displace the traps, the change in the condensate orbitals has a strong influence in the 
tunneling dynamics. Compared to the frozen model, which neglects this change, the full model 
gives a smaller frequency and higher amplitude of oscillation. We also found that in this latter 
case the full-model is much closer to its linear approximation than to the frozen model. 

As a final remark, it is clear that, by changing the initial condition, we can get a large variety 
of dynamical behaviors for the system. Our aim in this paper was to establish a framework 
where all these questions could be investigated in a scheme numerically feasible and which 
can be easily extended to a three dimensional calculation. The application of the formalism 
presented in this paper to this more realistic case is currently under investigation. 
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Figure Caption 

Figure 1. Plot showing time evolution of relative population (la), translation (lb-lc) and 
breathing (ld-le). The solid curve corresponds to the full model, the doted curve to the 
"frozen" model and dashed to the linear approximation. The initial condition is the relative 
phase at #(0) = 0.77T and the other variables are at its equilibrium values. The displacement of 
the trap zq is set equal to 0. The time is expressed in unit of ^ and the lengths in unit of z s h , 
which is the size parameter of the trap. See text for more details. 

Figure 2. The conventions are the same as in figure 1. The displacement of the trap z is set 
equal to 0.15z s h o - The initial phase is 0(0) = 0.77T and the other variables are at its equilibrium 
values. See text for more details. 
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Figure 1 
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